Modeling the Auxetic Transition for Carbon Nanotube Sheets 
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A simple model is developed to predict the complex mechanical properties of carbon nanotube 
sheets (buckypaper) [Hall et al, Science 320 504 (2008)]. Fabricated using a similar method to that 
deployed for making writing paper, these buckypapers can have in-plane Poisson's ratios changed 
from positive to negative, becoming auxetic, as multiwalled carbon nanotubes are increasingly mixed 
with single-walled carbon nanotubes. Essential structural features of the buckypapers are incorpo- 
rated into the model: isotropic in-plane mechanical properties, nanotubes preferentially oriented in 
the sheet plane, and freedom to undergo stress-induced elongation by both angle and length changes. 
The expressions derived for the Poisson's ratios enabled quantitative prediction of both observed 
properties and remarkable new properties obtainable by structural modification. 



PACS numbers: 62.25.-g, 61.46.Fg, 62.20.dj 

I. INTRODUCTION 

The ratio of percent lateral contraction to percent ap- 
plied tensile elongation is the Poisson's ratio. If the lat- 
eral dimension expands during stretching, the Poisson's 
ratio is negative and the material is called auxetic pQ. Re- 
cent interest in this counter-intuitive behavior originated 
from pioneering discoveries that partially collapsed foams 
and honeycombs [TJ [5], fibrillar polymers |3], and poly- 
mer composites [3] can be auxetic. This unusual property 
of auxetic materials results in various useful effects, such 
as increased indentation resistance and increased shear 
stiffness [5]. Possible or deployed applications of aux- 
etic materials are, for example, anti-ballistic vests, air 
filters, strain sensors, molecular-scale amplifiers, vascu- 
lar implants, gaskets, sound absorbers, artificial muscles, 
and wrestling mats [U 5, 6J. Due to their unusual and 
interesting properties, auxetic materials have been the 
subject of intense experimental and theoretical research 

[£113 E]. 

Recently, we showed that the Poisson's ratio of carbon 
nanotube sheets (buckypaper) can change from positive 
to negative as multiwalled carbon nanotubes (MWNTs) 
are increasingly mixed with single-walled carbon nan- 
otubes (SWNTs) [5]. While the in-plane Poisson's ratio 
for SWNT nanotube sheets was positive (about 0.06) and 
slightly changes until MWNT content reached 73 weight 
percent (wt. %), further addition of MWNTs decreased 
Poisson's ratio to —0.20 (Fig. 1). On the other hand, 
large positive Poisson's ratios were observed for the thick- 
ness direction: 0.33±0.14 and 0.75±0.30 for SWNT and 
MWNT sheets, respectively. A non-linear dependence of 
Young's modulus, strength, and toughness on MWNT 
content was also observed, though electronic conduc- 
tivity and density depended approximately linearly on 



MWNT content. A model incorporating an idealization 
of the complex structure of the SWNT/MWNT bucky- 
papers and the main deformation mechanisms was pro- 
posed in order to understand the behavior of the in-plane 
and thickness-direction Poisson's ratios within these nan- 
otube sheets [S]. We herein provide a full account of 
calculation methods and results, which were previously 
briefly outlined [9] . 

The Poisson's ratios of individual SWNTs and 
MWNTs have been theoretically obtained using analyti- 
cal [TO! HI [H] and atomistic models including empirical 
potentials [T3J E], tight-binding-based approaches [TS] . 
and ab initio methods [16] . All of these investigations 
predict positive values for the Poisson's ratios. Calcula- 
tions using density functional theory [16 provide values 
between 0.12 and 0.19, comparable to the Poisson's ra- 
tio for the basal plane in graphite (0.16) [T71 HH]. How- 
ever, when individual nanotubes are assembled together 
in sheets containing fiber networks, either negative or 
positive Poisson's ratios can arise [5|. Analyzing carbon 
nanotube sheets, Berhan et al. [19 showed using Euler 
beam-network simulations that increasing the number of 
interfiber connections can lead to improvements in car- 
bon nanotube sheet stiffness. 

Fiber networks have been subject of extensive inves- 
tigation over the last six decades (e.g. ^ZU\ [5TJ [551 H31 
l24l [25]). Using effective-medium theory, Cox's pioneer 
work on fibers [20] predicts the effective moduli of two- 
dimensional fibers considering only the fiber axial defor- 
mation. Analyzing a solid mat of fibers, Cox predicted 
the possibility of a negative Poisson's ratio in plane of 
the paper existing concurrently with a high positive value 
through the thickness of the paper. However, the Pois- 
son's ratio measured in-plane for ordinary fiber networks, 
like writing paper, are rather large and positive. Exten- 
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sions of the Cox's model have been developed where fiber 
bending, elongation, and contraction, as well as Poisso- 
nian distribution of fiber segment lengths, are taken into 
account [2TJ [53J H3] . Producing carbon nanotube net- 
works from dispersed aqueous nanotube suspensions [3], 
we were able to provide an experimental realization of 
Cox's prediction that a negative Poisson's ratio can exist 
for paper-like fiber mats. The present goal is to describe 
a simple model that provides a realistic, though simpli- 
fied, description of the structural nature of buckypaper, 
and use this model to explain why the observed Poisson's 
ratio of multiwalled and single walled carbon nanotube 
buckpaper differ in sign. 

This paper is divided into the following: Section II 
briefly describes the experimental details of the buckypa- 
per fabrication and the methods used to determine the 
Poisson's ratios. The model and the derivation of the 
Poisson's ratios are presented in Sec. III. In Sec. IV, we 
present the results obtained from the model and discuss 
the main results in Sec. V. Conclusions and final remarks 
of this work are summarized in Sec. VI. 
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acid-reflux purified SWNTs (Buckypearls) differed in in- 
plane Poisson's ratio by only 0.034. The nanotube bucky- 
paper was fabricated by vacuum filtration of an ultrason- 
ically dispersed aqueous nanotube suspension containing 
Triton X-100 surfactant, washed with successively wa- 
ter and methanol, vacuum drying (85°C for two days), 
and then peeling the nanotube sheet from the filter [3"U] . 
Nanotube sheet thickness was held approximately con- 
stant (50 fim). Scanning electron microscopy indicated 
that SWNTs and MWNTs are intimately commingled 
in sheets comprising both nanotube types. The aver- 
age angle between the nanotube length direction and the 
sheet plane (41.7° for MWNT sheets and 45.0° for SWNT 
sheets) was determined by diffraction for incident x-rays 
in an in-plane sheet direction, using the dependence of 
diffraction intensities on azimuthal angle [31] ■ Fig. 2 
shows scanning electron microscopy images of buckypa- 
per surfaces for sheets produced with different MWNT 
wt. % contents. 

Reported in-plane mechanical properties measure- 
ments are for 2 mm x 12 mm carbon nanotube strips 
during deformation at 0.10% strain/minute. Poisson's 
ratio measurements utilized nanotube sheets coated with 
trace Ti02 particles for marking position. Digital images 
were captured during constant rate tensile deformation, 
and interpreted using image correlation software (Vic-2D 
Correlated Solutions, Inc., West Columbia SC, USA) to 
obtain changes in the separations between thousands of 
Ti0 2 particles as a function of tensile stress, correspond- 
ing sheet strains in stretch and lateral directions, and the 
Poisson's ratio. The thickness-direction Poisson's ratio 
was obtained from scanning electron micrographs show- 
ing sheet thickness versus applied in-plane tensile strain. 
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FIG. 1: Measured in-plane Poisson's ratio vs. MWNT content 
in SWNT/MWNT sheets. The continuous line is only a guide 
to the eye. 

II. EXPERIMENTAL DETAILS 

The utilized MWNTs (-12 nm in diameter, -200 ^m 
long, and containing about 9 walls) were harvested from 
nanotube forests prepared by chemical vapor deposition 
(CVD) of acetylene gas [26]. The SWNTs (-1.0 nm di- 
ameter and below micron long) j27j [28] were commer- 
cially produced by Carbon Nanotechnologies, Inc. using 
the HiPco synthesis method by high pressure CVD of 
carbon monoxide gS]. The MWNTs had below 2 wt. % 
catalyst while the unpurificd HiPco nanotubes have high 
wt. % catalyst, though low volume % catalyst. Cata- 
lyst concentration is apparently unimportant for the used 
as-synthesized HiPco SWNTs, since buckypaper sheets 
made from commercially obtained as-synthesized and 



III. DEVELOPMENT OF THE MODEL 
Structural model 

As we can see from Fig. 2, the structures of the carbon 
nanotube sheets are very complex, nanotubes and nan- 
otube bundles having various diameters meander in three 
dimensions, like cooked spaghetti. Therefore, the chal- 
lenge is to formulate a model for these intractably com- 
plex nanotube sheets that is sufficiently simple that it can 
be tested and used for prediction of future results. A first 
approach to achieve this is to simplify the complex mor- 
phology of the meandering nanotubes in the thickness 
direction to a zigzag set of struts as illustrated in Fig. 3 
(a). The bends at the zigs and zags enable tractable rep- 
resentation of the observed deviation of nanotubes from 
perfect in-plane alignment and the geometrical effect of 
nanotube straightening on in-plane and thickness direc- 
tion Poisson's ratios. The bend force constant at the zigs 
and zags correspond to the effective force constant for 



FIG. 2: Scanning electronic microscopy images of the surfaces 
of buckypaper containing (a) wt. %, (b) 47.1 wt. %, (c) 
72.7 wt. %, and (d) 100 wt. % MWNT content. Different 
magnifications are shown in top and bottom parts of each 
image. The scale bars for the upper and lower images in (a)- 
(d) correspond to 2^tm and 200 nm, respectively. 



FIG. 4: (Color online) Schematic view of an "egg rack"-type 
structure representing the real complex structure of the nan- 
otube sheets. The arrow indicates the position where two 
model fibers would be in contact. Due to the connectivity 
of the structure, tensile loading in the in-plane direction (xy 
plane) produces a negative Poisson's ratio in the plane of the 
structure and a positive Poisson's ratio in the thickness direc- 
tion (xz and yz planes) if the only deformation made is angle 
bending. 



elongating a meandering nanotube in a network of inter- 
acting neighboring nanotubes. Coupling between inter- 
secting nanotubes is at junctions, where the zags from 
one layer of zigzag chains are coupled to the zigs for the 
next layer of chains. 




FIG. 3: (Color online) Schematic representation of the ap- 
proximation used to derive a simple model for the nanotube 
sheets, (a) A SWNT bundle or a MWNT that is deformed 
along the in-plane and thickness directions (top) is repre- 
sented by a set of struts forming a zigzag arrangement (bot- 
tom), (b) The complex nanotube sheet morphology (top) 
is represented by a simplified model of an ordered structure 
(bottom). 



The second simplification is to describe the disordered 
observed structure by an mechanically equivalent struc- 
ture comprising layers that are periodic within the sheet 
plane as an ordered structure (Fig. 3 (b)). A simple 
example of an ordered structure that could be used for 
representing the nanotube sheets is similar to an egg- 
rack. This structure is composed of oppositely facing 
'four-legged claws' arranged on a square grid. Fig. 4 
depicts a schematic view of a "egg-rack" -type structure. 
Grima et al. |35j demonstrated that when this structure 
is loaded in tension, the connectivity of the claws forces 
them to open in all directions, hence producing a nega- 
tive Poisson's ratio in the plane of the structure and a 
positive Poisson's ratio in the thickness direction. These 
two features of such type of structure are observed in 
the real MWNT sheets. The meandering of fibers can be 
represented in this structural model as indicated in Fig. 
4 by the magenta and red struts. Moreover, this repre- 
sentation also permits the fibers to cross each other, like 
they do in the real structure (Fig. 2). 

While the egg-rack model incorporates aspects ob- 
served in real nanotube sheets, such as preferential nan- 
otube orientation in the sheet plane (but with positive 
and negative deviations from in-plane orientation) and 
quite different Poisson's ratios for sheet plane and sheet 
thickness directions, it lacks the needed isotropy for in- 
plane mechanical properties. Anisotropic models have 
been previously deployed to represent sheets that have 
isotropic in- plane properties [36]. Uncertainty in the sign 
of Poisson's ratio resulted from the need to average in- 
plane properties to obtain predicted properties for sheets 
having isotropic in-plane mechanical properties. Thus, in 




FIG. 5: (Color online) (a)-(b) Three-dimensional structure 
used to represent nanotube sheets. Each strut represents a 
fiber and each ball a fiber junction. In (b) the <f) angle is 
represented through virtual struts (yellow dashed lines), (c)- 
(d) Lateral views of the structure with different inclination 
angles (7) of the struts with respect to the sheet plane: (c) 
20° and (d) 40°. 

order to avoid the averaging process, a structural model 
showing in-plane isotropy is desired. 



The simplest model that provides all these key features 
of the nanotube sheets has the hexagonal space group 
P6222 and inter-nanotube non-covalent junctions located 
at (0.5, 0, 0) and equivalent locations in the unit cell. In 
this model, the nanotubes (and nanotube bundles) are 
represented by zigzag chains parallel to the sheet plane 
(with angle between the struts and the basal plane of ±7 
and an inter-strut angle of 29 = it — 2j) . Zigzag chains 
in one nanotube sheet layer connect non-covalently with 
those in the next layers at the extremes of the zigs and 
zags, where torsion about the contact enables change in 
the intersection angle <j) between nanotubes. A three- 
dimensional view of the structure is presented in Fig. 5. 
The same mechanical properties result for the closely re- 
lated structure shown in Fig. 6, in which each successive 
layer of zigzag chains are equally likely to be added in 
either of two possible directions. Before deformation by 
application of tensile stress along the ^-direction, struts 
1 have length L 10 and length projected onto the sheet 
plane (xy plane) of Li p o, and struts 2 have the same 
length (L20) and projected length (L 2p o) as for struts 1 
(Fig. 6). For the special case of 7 = 0, Li pQ = L w . Thus 
the structural model depicted in Fig. 6 is the approxi- 
mation used here to represent the complex structure of 
nanotube sheets. 

The last step on the model definition is to decide which 
nanotube deformation mechanisms in the real nanotube 
sheets will be taken into account within the model. Due 
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FIG. 6: (Color online) (a)-(b) Labeled schematic illustration 
of the hexagonal model structure viewed perpendicular to the 
sheet plane, where the zigzag nanofibers in the same layer 
share the same color, (c) Labeled lateral view of a zigzag 
nanofiber in this structural model. 



to the complex morphology of the nanotube sheets, the 
inclusion of all possible deformation types would not be 
possible for a predictive model. Therefore, we incorpo- 
rated only the deformation modes we considered to be 
most important for describing the mechanical properties 
of the nanotube sheets. Thus, the following elastic defor- 
mations were considered: nanotube (or nanotube bun- 
dle) axial stretching, represented by the force constant 
k s ; nanotube bending due to changes in the 9 angle, rep- 
resented by the bending force constant kg; and due to 
changes in the torsional angle <fi between coupled inter- 
secting nanotubes, represented by the torsional force con- 
stant kt- These force constants are effective values, aris- 
ing in the complex real structure from the energy needed 
to straighten meandering nanotubes and change the an- 
gle between intersecting nanotubes. 

Before we proceed with the determination of the ex- 
pressions for the Poisson's ratios of the hexagonal model, 
it is useful to derive the effective force constant fc s fc of 
a single zigzag nanofiber when both strut stretch and 9 
angle deformation are included. The zigzag chain is de- 
picted in Fig. 7. The dimensions of the unit cell are 
given by L x = 2Lsn\9 and L z — Lcos9, where L is the 
length of the strut and 29 is the inter-strut bond angle. 
The zigzag chain force constant is related to the applied 
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force F through k sb = 2FjAL, 




FIG. 7: (Color online) Zigzag chain (top) and its correspond- 
ing unit cell (bottom). 

The force F can be written in terms of the parallel F\ 
(= Fsin6>) and perpendicular F 2 (— Fcos9) components 
with respect to the strut. The forces F\ and F 2 will then 
produce the corresponding displacements 8\ = F\/k s and 
82 = Fz/kg. Thus, the variation AL X can be obtained 
using the components of Si and <5 2 , leading to AL X = 
2F(sin 2 8/k s +cos 2 8/kg). Consequently, k sb will be given 
by: 

k sb = . (1) 

kg sin 6 + k s cos 2 9 

We can see that when kg/k s — > 0, k sb — ► kg/ cos 2 and 
when kg/k s — * 00, k sb — > fc s /sin 2 #. 

Derivation of the Poisson's ratios 



constraint that all layers have the same tensile-direction 
and width-direction strains: (i) 2 + 20i = ^ and (ii) 
Li p = 2L 2p sin(0 2 /2). Due to constraint (ii) ALi p is 
written as 

ALi p = L cos((0 o + A0 2 )/2) A0 2 + 2 sin(0 o /2) AL 2p . 

(3) 

where we used L\ p q = L 2p o = Lq. In the elastic regime 
A0 2 is small, therefore, cos((0o + A0 2 )/2) — v3/2 — 
A0 2 /4, thus ([3]) can be written as 

where <5 = AL 2p /L - Since we are comparing terms up 
to the first order in 5 and e in the above equation, it 
is reasonable to assume that (A0 2 ) 2 is very small when 
compared to e and S. Therefore, the constraints (i) and 
(ii) are expressed as: (a) A0 2 + 2A0i = and (b) e = 
y/3A<f> 2 /2 + 5, with e,8 < 1. Using (b) the energy E 
becomes 

E = (k sb /6)[(eL ) 2 + 2(SL ) 2 ] + (3/2)(fc 4 /6)(L o A0 2 ) 2 . 

(5) 

Defining / = QE / (ktLify and R = k sb /kt the function 
to minimize can be expressed as 

/-(<5,A0 2 ) = i?(e 2 + 2£ 2 ) + (3/2)A0 2 , (6) 

subject now to the following constraint: 

5 (<5,A0 2 ) = e-^A0 2 -5. (7) 

The solutions can be obtained by using the Lagrange's 
multipliers method, solving the system: 



When an in-plane tensile stretch is applied along the 
x-direction, causing a small tensile strain e = AL\ p / L\ p q 
(<S 1), strut lengths and angles change to L>i p — Li p0 + 
ALi p , L 2p = L 2p o + AL 2p , 0i = 0io + A0i, 2 = 
2O + A(f> 2 (0io = 020 = 0o = 60°), and likewise for 
8. The total energy per strut (E) needed for producing 
a given small tensile strain in terms of angle bend and 
strut length changes is given by: 



E = (fc s& /6)[(ALi p ) 2 + 2(AL 2p ) 2 ] 
+ (fc t /6)[2(Li pO A0i) 2 + ( J L 2pO A0 2 ) 2 ] 



(2) 



Minimization of E provides all changes in lengths and 
angles for a specified small tensile strain in the in-plane 
direction. From these changes, the in-plane (yi) and 
sheet thickness direction (^3) Poisson's ratios can be ob- 
tained. The energy must be minimized subject to the 



(xaa,\ n d/(<5,A0 2 ) gMA0g) 
3(d,A0 2 ) = O, — =0, — =0,(8) 



9A0 2 



where /(<5,A0 2 ) = f(8,A<j) 2 ) - Ag(<5,A0 2 ) and A is the 
Lagrange's multiplier. The solution for this system is 



2R 

A0 2 = -= e 

V3(l + R) 



and 



(9) 



(10) 



1 + R 

The energy at the minimum will then assume the value 



E,, 



R(R + 3) , t2 2 
6(1 + ii) * ° £ 



(11) 
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The in-plane Poisson's ratio v\ is given by 
/ AL*/L* 



\ALi p /L 



(12) 



where L* = L>2 P cos (02/2). Neglecting second order 
terms we can write 



AL* ~ sin(0 o /2)A0 2 + cos(<fo/2) AL 2p , (13) 



or 



Consequently: 



3(1 + i?)' 



(14) 



(15) 



Following the definitions of Fig. 7, the Poisson's 
ratio for a single zigzag chain is given by v Z i g — 
— (AL Z /L Z )/(AL X /L X ). The term AL Z is calculated us- 
ing Si and (>2, leading to AL Z = F sin 8 cos 8(1 / k s — l/ kg) . 
Therefore, 



tan 2 8(k s /k e - 1) 
tan 2 8 + k a /he 



(16) 



The thickness-direction Poisson's ratio vz is deter- 
mined through a three-layer average, i.e., 



^ = (1/3) - 



(i)/ r (i) 



AL Z ' /L 
ALip/Lo 



A/ -(2) /r (2)\ 

< (17) 



where AL Z 1 ^ / refers to the thickness change of the 
strut i due to change in the tensile direction change 
AL\ v j Lq- 

Manipulation of (17) leads to: 



vz = 



1 ( alP/l™ 

AL\ v j Lq 



AL Z 2) /L Z 2) AL 2p jU 
AL<2 P /Lq ALip/Lo 



(18) 



or 



2u 7 _ 



AL-2 P /Lo 
3 \'*' y ' ~"^ g AL lp /L 

which finally yields to 



^v zig I 1 + 2 



3 + R 
3(1 + R)' 



(19) 



(20) 



Using the following definitions for the ratios between 
force constants s = k t /k s and r = kg/kt, the Poisson's 
ratios can be expressed as 



3 + (3 



and 



(1- S )(1 + P) 
(tan 2 7 + s)(3 + /3)' 



(21) 



(22) 



where (3 = 3k t /k s b = 3[1 + (s — 1) cos 2 j]/r. 

For the case where the nanotube struts have effectively 
infinite modulus (,s = 0), the expressions reduce to 



1-/3 
3 + /3 



and 



vz 



1 + /3 j 
3+^3 ° 0t 7 ' 



(23) 



(24) 



but now with (3 — (3sin 2 7)/r. Using a different model, 
which includes a host of structural and force constant 
parameters in (3, the above dependence of V\ on (3 has 
been predicted for sheets of cellulose-based papers 38J. 



IV. MODEL RESULTS 



Poisson's ratio behavior 



We can see from expressions (211 and (22 1 that when 



changes in the angle between intersecting nanotubes are 
negligible in comparison to changes in nanotube length 
(due to both nanotube stretching and changes in the 
angle 8), expressed as (3 — > 00, the most negative in- 
plane Poisson's ratio (v\ — > —1) and most positive thick- 
ness direction Poisson's ratio (vz — > cot 2 7) are ob- 
tained. On the other hand, when changes in the an- 
gle between intersecting nanotubes are much easier than 
changes in the nanotube length ((3 — » 0), v\ — > 1/3 and 
vz — > (1/3) cot 2 7. We can see that in both limits the 
thickness-direction Poisson's ratio is positive while the 
in-plane Poisson's ratio undergoes a transition from neg- 
ative (auxetic) to positive (non-auxetic) values. 

This non-auxetic /auxetic behavior can be simply visu- 
alized by noting that two neighboring nanotube layers in 
the Fig. 6(a) are coupled like the struts of a wine rack. If 
rotation between struts dominates, like for an ordinary 
wine rack, the Poisson's ratio is positive. If this tor- 
sional rotation of struts is blocked (by welding together 
the struts) and the struts are stretchable but not bend- 
able, strut length increases produce a negative Poisson's 
ratio. 

These results indicate that the present model allows a 
qualitative description of the experimental behavior ob- 
served for the carbon nanotube sheets [5] if the mixing 
of MWNTs in SWNT sheets can be represented by an 
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effective change on the ratios r and s as MWNTs are in- 
corporated, fn order to estimate the utility of the model 
we can determine the r and 7 values (considering the 
approximation s = 0) that provide the observed values 
of Poisson's ratio (by matching the expressions ( p3| and 
(24 1 with the experimentally observed Poisson's ratios). 
The expressions ( 23 1 and ( 24 1 can be inverted yielding to 



3(1 -v\) 



(3^-1x^-1-2^ 



and 



cos 7 



2^3 



1 — V\ + 2i/ 3 



(25) 



(26) 



Using the experimental data for MWNT sheets {y\ — 
—0.20 and v 3 = 0.75) the derived values from the above 
expressions are r ~ 0.67 and 7 ~ 42°. Similarly, for 
SWNT sheets (v x = 0.06 and v z = 0.33) we have r ~ 2.28 
and 7 ~ 50° . The predicted 7 values are consistent with 
average angles from x-ray diffraction of 41.7° for MWNT 
sheets and 45.0° for SWNT sheets [9]. With these values 
of r and 7 we have that the ratio of the (3 parameters for 
MWNT and SWNT sheets is /3 M wnt//%wnt - 2.6. 



Negative linear compressibility 



alignment might be obtainable by either using high pres- 
sure for the filtration step used for sheet fabrication or by 
annealing the as-fabricated sheets under high mechanical 
load. 
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FIG. 8: (Color online) The relationships between in-plane 
and thickness-direction Poisson's ratios as a function of indi- 
cated average nanotube orientation angle 7 for SWNT and 
MWNT sheets having the force constant ratios r that yield 
the measured Poisson's ratios (circles). The straight (blue) 
line represents the equation 1 — v\ — vz = 0. Values of 7 on 
the right side of this line are predicted to produce negative 
linear compressibilities for SWNT and MWNT sheets. 



Negative Poisson's ratios are sometimes accompanied 
by much rarer mechanical properties: negative lin- 
ear compressibilities and negative area compressibility, 
meaning that a material expands in either one or two or- 
thogonal directions when hydrostatic pressure is applied 
|37j . A negative linear compressibility is the inverse of 
another strange property - increasing density when elon- 
gated in a direction where linear compressibility is neg- 
ative, and both require that 1 — v\ — v 3 < |37j . Using 



( 25 1 and ( 26 ) we have that this condition becomes 



1 - V\ - 2/3 = 



2r — 3 cos 2 7 — r cot 2 7 + 6 sin 2 7 
3(r + sin 2 7) 



(27) 



Thus negative in-plane compressibility, negative area 
compressibility for the sheet plane, and stretch densi- 
fication is predicted for cos 7 > y/2/3, which implies 
7 < 35.3°. Fig. 8 shows the predicted behavior of v\ 
against ia, for the pure SWNT (r = 2.28) and MWNT 
(r = 0.67) samples for the case where s — 0. The average 
7 needed for achieving these properties will decrease as 
a result of in-plane nanofiber meandering, since only the 
tensile strain component resulting in thickness change ef- 
fects V3. Since the enhanced degree of in-plane alignment 
needed to realize negative linear compressibilities is not 
large, the needed improvement on the degree of in-plane 



V. DISCUSSIONS 

In order to discuss the role of the nanoscale con- 
stituents of nanotube sheets on the behavior of the in- 
plane Poisson's ratios we will proceed as follows. Firstly, 
we will show that beam bending in response to tensile 
stress presents the same effect for increasing Poisson's 
ratio as torsional rotation. This allows the use of the 
force constant ratio between strut bending to stretching 
for understanding the in-plane Poisson's ratio values for 
nanotube sheets. Secondly, we will estimate this ratio for 
the nanotube sheets based on experimentally accessible 
quantities and use this result for predicting the in-plane 
Poisson's ratios for SWNT and MWNT buckypaper. 



Nanotube bending vs. nanotube stretching 

In analogy with a wine rack, positive in-plane Pois- 
son's ratios would result if (i) the hinges are welded to 
struts to prohibit torsional rotation and (ii) the struts 
are much easier to bend than to stretch. Nanotube 
beam bending in response to a tensile stress within the 
sheet plane changes the effective angle between intersect- 
ing nanotubes, and produces a corresponding increase in 
Poisson's ratio, similarly to the response when there are 
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changes in torsional angle for the model of Fig. 6(a). If 
fiber beam bending is the predominant deformation that 
changes the effective angle between intersecting fibers 
and fiber deviation from in-plane orientation is neglected 
(7 = 0), (3 = 3kb/k s b, where the force constants for fiber 
bending and tensile fiber elongation are kb and k s b, re- 
spectively. 

In order to obtain the dependence of v\ on kb/k s b we 
proceed as follows. Since only the in-plane Poisson's 
ratio v\ is presently being evaluated, we used a two- 
dimensional (2D) sheet structure to represent the depen- 
dence of v\ on effective kb/k s b- These 2D sheets look like 
the projection of the structure shown in Fig. 6(a), ex- 
cept that six co-planar struts meet in 2D at each junction. 
With the purpose of excluding torsional angle changes, 
but enabling strut length changes and strut bending, the 
struts were represented by long chains composed of N 
"atoms" (N large), which were allowed to undergo bond 
angle bending at each atom and elongation of the bonds 
between atoms (see Fig. 9). The struts were connected 
to artificially contracted six membered rings that were 
so small and so rigid with respect to angle and dimen- 
sional changes that they acted as an junction that does 
not allow the equivalent of torsional rotation. For a sin- 
gle value of kb/k s b the structure was geometrically op- 
timized and the Poisson's ratio was obtained from the 
stress tensor derived from the second derivative of the en- 
ergy with respect to the strain. These calculations were 
carried out using the Cerius 2 open force field molecular 
mechanics [j^S HD] ■ By varying strut bending and strut 
elongation moduli arbitrarily, the calculations provided 
v\ = (1 — /3)/(3 + /?), with [3 — 3kb/k s b as shown in Fig. 
10 [40] . An illustration of the deformations presented by 
the 2D model under tensile strain is depicted in Fig. 11. 
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FIG. 10: (Color online) In-plane Poisson's ratio vs. kb/k a b 
obtained from molecular mechanics calculations. 



Nanoscale origin of the auxetic property of 
nanotube sheets 

Having shown that beam bending has a similar effect 
for increasing in-plane Poisson's ratio as torsional rota- 
tion, indicating that the ratio kb/k s b plays an important 
role in the understanding of the in-plane Poisson's ra- 
tio behavior, we will estimate kb/k s b for the nanotube 
sheets. Before that, it is worth mentioning that when 
7 = [9 = 90°) and changes in the 9 angle are negligible 
compared to changes in the strut lengths (kg/k s — ► 00) 
k s b — > k s . Thus, the differences between MWNTs and 
SWNT bundles will be associated with the intrinsic me- 
chanical properties of each of them, being represented 




Strut length 



FIG. 9: (Color online) Schematic representation of the 2D 
sheets used in molecular mechanics calculations, (a) Repre- 
sentation of the unit cell, (b) Detailed view of a region of (a) 
showing the main spring constants used to describe the 2D 
deformations |4U] . 




FIG. 11: (Color online) Resulting deformations of the 2D 
sheet model under tensile strain for (a) k s /k s b — > (non- 
auxetic) and (b) k s /k s b — > 00 (auxetic). Top (bottom) dia- 
grams correspond to the underformed (deformed) state. 
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here by effective spring constants related to nanotube ax- 
ial stretching k s and to the nanotube lateral bending kb- 
Therefore, the ratio kb/k s will be different for each nan- 
otube type and will depend upon their nanoscale struc- 
ture. 

Consider a perfect, straight SWNT of radius r t hav- 
ing length I described by the elastic tube model [T31I34"] . 
In this case, k s = Y s A/l and kf, = 3YbI/l 3 , where Y s 
and Y b are the Young's modulus for stretching and bend- 
ing, respectively, A is its cross-section area, and / is 
the moment of inertia. Using A = 2irr t h, Y s = C/h 
and Yf,I = K = TrCr 3 [33], where h is the graphene 
sheet thickness, we obtain k\,jk s — (3/2)(r t /7) 2 . Here we 
assumed that bending and stretching Young's modulus 
have the same value. 

For a bundle of radius Rb composed by p SWNTs of 
radius r t that can freely slip along each other, acting 
independently, the bending force constant is expressed 
as pkt,. Assuming that the load is carried by the SWNTs 
on the perimeter of each bundle |H] , the stretching force 
constant will be given by k s = (C/h)(A cxt /l), where A ext 
is the cross section area of the outer nanotubes in the 
SWNT bundle. For large diameter bundles (Rb ^> r t ) 
A C xt — 2irRb(2r t ). For SWNTs closely packed into a two- 
dimensional hexagonal array, the volume fraction V p is 
0.906 thus p = Vb{TiRl)/{^) = V b (R b /r t ) 2 . Therefore, 



bundle 



3hV b f R b 

l 2 



(28) 



The bending stiffness of a MWNT can be expressed 
simply as a sum of bending stiffness of the individual 
and independent SWNTs [33], i.e., K = 7r£"i n . C i r \j\> 
where i is the corresponding wall number of the MWNT. 
For a large diameter MWNT it is reasonable to assume 
that all Cj for this MWNT are the same and equal to 
the in-plane stiffness C. Considering that the MWNT 
is composed only by armchair SWNTs then r t u-\ = ai, 
a = 0.68 A thus K = <Ca 3 , where £ = * 3 - 
Assuming that the load is supported by only the outer 
wall of the MWNT then k s = (C/h)(A outer /l) where 
Amtcr = 2irr t ( nf )h. Consequently, 



MWNT 



3a 3 

~2~ 



r- 



r t (n f ) 



(29) 



From eqs. (28 1 and (29) we can see that to de 



crease the in-plane Poisson's ratio (through the increase 
(3 = 3k b /k s ), inter-junction lengths / should decrease 
or, equivalently, sheet density should increase. For 
SWNT sheets, increasing bundle radius also leads to a 
decrease of the Poisson's ratio. For the limiting case 
where the bundle is formed by only one SWNT, kb/k s — 
(3/2)(r t /l) 2 . The largest geometrically possible value of 
r t /l is sin(60°)/2, which corresponds to the physically 
unreasonable case where each layer within the nanotube 



sheet comprises straight nanotubes that are close packed 
within the layer. From (23) and f3 = 3kb/k s , this hypo- 
thetical buckypaper of perfectly straight, infinitely long, 
unbundled SWNTs can not have a Poisson's ratio be- 
low 0.04. Since for buckypaper-like sheets comprising 
long circular solid fiber is the same as for a SWNT when 
the effective Young's modulus for bending equals that 
for tension, the predicted v\ is also 0.04 or higher, and 
likely much higher since buckypaper-like sheets or ordi- 
nary paper do not have fibers that are close packed in a 
plane. 

For MWNT sheets, a decrease of the in-plane Poisson's 
ratio can be obtained by increasing the number of inte- 
rior walls (which increases £) and the MWNT diameter. 
While all nanotube walls contribute additively to kb, only 
the outer wall contributes to k s unless the MWNTs are 
extremely long. However, the effects of these structure 
changes are not simple, since increasing kb and decreasing 
I can decrease nanotube meandering between junctions 
and this decrease of meandering can provide a positive 
contribution to k s . 

In order to have a rough estimate of kb/k s for MWNTs 
and SWNT bundles present in the nanotube sheets (Fig. 
2), we use experimental data [9] showing that the SWNTs 
have an average diameter of about 1.0 nm and a average 
bundle diameter of 20 nm and that the MWNTs have an 
outer diameter of about 12 nm, and contain about nine 
walls (Fig. 12). From these data rii — 48, rif = 88, 
£ = 3135888, and k b /k s ~ 2/1 2 and k b /k s ~ 247/Z 2 
(I in nm) for SWNT and MWNT sheets, respectively. 
For MWNTs and SWNT bundles sharing approximately 
the same length we can see that /3mwnt//5swnt ~ 100, 
two orders of magnitude larger than the previously ob- 
tained value of 2.6 when the nanotubes are considered 
completely rigid. While these estimates lead to positive 
and negative values for the in-plane Poisson's ratio for 
SWNT and MWNT sheets, respectively, in qualitatively 
agreement with experiment, quantitative agreement is 
poor (0.3 and -0.6 for SWNT and MWNT sheets, re- 
spectively, with I ~ 10 nm). As we will see, this differ- 
ence lies on the assumption that the effective force con- 
stant for elongating the nanotubes corresponds to the 
modulus of an individual straight nanotube. Because of 
nanotube meandering, and possibly elasticity at inter- 
nanotube junctions, this assumption is not valid. 

In order to better estimate v\ we can proceed as fol- 
lows. The Young's modulus of the nanotube sheet is 
obtained by taking the second derivative of the per-strut 
energy of (11) with respect to e, 



Y = 



L k sb (l - ui) 
2V S 



(30) 



where we used R = k s b/k t = 3(1 + ^i)/(l — 3v\) from 
(15) and V s represents the volume per strut in the sheet 
structure. With this expression and the experimentally 
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observed Young's modulus for the nanotube sheets, we 
can obtain a more precise estimate of k sbl instead of only 
k s , as previously considered in the limiting case of 7 = 0. 

Deriving the effective strut lengths for the highly dis- 
ordered SWNT and MWNT buckypaper is important, 
since Lq and V s are needed for comparing theory with 
experiment. In order to calculate these parameters, we 
look at the intersecting nanotubes (or nanotube bundles) 
as being stacked in the thickness direction like layers of 
logs having an effective diameter D, where D is the sum 
of the covalent diameter of the nanofiber and the 0.34 nm 
van der Waals diameter of carbon. Correspondingly, the 
volume per strut is DL\ sin(120°). If the strut weight 
per strut length is Wl and the measured nanotube sheet 
density is p, then p = W L /[DLl sin(120°)]. Using the 
observed densities for the MWNT (0.343 g/cm 3 ) and 
SWNT sheets (0.692 g/cm 3 ) jS], the corresponding calcu- 
lated L for MWNT and SWNT struts are 54.3 nm and 
39.5 nm, respectively. While these distances seem shorter 
than suggested by the micrographs of Fig. 2, note that 
these micrographs are for the sheet surface (the face orig- 
inally in contact with the filter membrane) and do not 
provide the junction density and corresponding Lq in the 
buckypaper interior. Using this strut volume, eq. (30) 
becomes 



k s b — 



2YD sin(120°) 



(31) 



The k b for the MWNTs is obtained by the sum of 
bending force constants for all component SWNTs (with 
C = 345 J/m 2 [33]) and the k b for SWNT bundles is 
derived from the measured average Young's modulus for 
bending (Y b ) 20 nm diameter SWNT bundles (50 GPa) 
[4"3"] , using the force constant for bending a solid cylin- 
drical rod k b = 3nrfY b / (4/ 3 ) Using the observed 
Young's modulus (1.81 GPa and 3.21 GPa for MWNT 
and SWNT sheets, respectively) (=k b /k sb ) can be self- 
consistently obtained from the relations above, which 
provides 1.84 and 0.42 for MWNT and SWNT sheets, 
respectively. Therefore, the predicted in-plane Poisson's 
ratios are —0.17 for MWNT buckypaper (vs. the ob- 
served -0.20) and 0.17 for SWNT buckypaper (vs. the 
observed 0.06). Considering Y b — 81 GPa, which is 
within the range of experimental uncertainty [43] . the 
calculated v\ for SWNT buckypaper assumes the ob- 
served value. For this latter estimate we can see that 
/3mwnt//?swnt — 2.4 in agreement with the previously 
predicted value of 2.6 for the case where spring constant 
for nanotube stretching is much larger than strut tor- 
sional rotations (s — k t /k s — ► 0), fitted from observed 
in-plane and thickness direction Poisson's ratios. 

Having predicted the Poisson's ratios for the MWNT 
content limiting cases of wt. % (SWNT sheets) and 100 
wt. % (MWNT sheets) we expect that MWNT/ SWNT 
mixing in buckypapers can be represented by altering the 
intrinsic characteristics of the fibers (nanotube bending 




12 nm 



20 nm 



FIG. 12: Schematic representation of the (a) SWNTs and (b) 
MWNTs present in nanotube sheets (Fig. 2). 



and elongation force constants) as well as by modifying 
the structural morphology (sheet density, fiber connec- 
tions) to provide the observed intermediate Poisson's ra- 
tios values. Further investigations are necessary to ex- 
plore all these possibilities and thus to provide a better 
understanding of the behavior of the Poisson's ratios for 
buckypapers having intermediated MWNT content. 



VI. SUMMARY AND CONCLUSIONS 

We proposed a relatively simple model to represent 
the complex morphology of carbon nanotube sheets and 
to describe the behavior of the Poisson's ratios. The 
model incorporates key structural features of the nan- 
otube sheets: isotropic in-plane mechanical properties, 
nanotubes preferentially oriented in the sheet plane, and 
freedom to undergo stress-induced elongation as a result 
of straightening meandering nanotubes and changing the 
angle between intersecting nanotubes. The nanoscale ori- 
gin of the constituent elements of the nanotube sheets is 
shown to play a fundamental role in determining the me- 
chanical behavior that leads to auxetic features, specially 
the ratio between the force constants associated with the 
fiber bending and elongation. Qualitative and quanti- 
tative agreement with experiment were obtained using 
the Poisson's ratio expressions derived from the model, 
encouraging its use in predicting future results and new 
properties. Specifically, negative linear compressibility 
(material expansion in the sheet plane under hydrostat- 
ically compression) was predicted for carbon nanotube 
sheets with the average angle of the fibers with respect 
to the sheet plane smaller than about 35° , a value which 
is not too far from already produced nanotube sheets. 
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